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If many micelles adsorb onto the same polymer molecule then they are said to form a necklace. A 
minimal model of such a necklace is proposed and shown to be almost equivalent to a 1-dimensional 
fluid with nearest-neighbour interactions. The thermodynamic functions of this fluid are obtained 
and then used to predict the change in the critical micellar concentration of the surfactant in the 
presence of the polymer. If the amount of polymer is not too large there are two critical micellar 
concentrations, one for micelles in necklaces and one for free micelles. 
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I. INTRODUCTION 



Experiment has shown that some micelles adsorb strongly onto polymer chains. If the polymer molecule is large 
enough many micelles can adsorb onto one molecule, forming a necklace of micelles along the polymer chain. The 
most widely studied example is the adsorption of micelles of sodium dodecyl sulphate (SDS) onto poly(ethylene oxide) 
(PEO), see Refs. and references therein. See Ref. H) for work on other systems. The interaction energy of a 
micelle with a polymer chain can be strong, many times kT the thermal energy, and has a large effect on both the 
polymer and the micelle. The favourable interaction of a micelle with the polymer reduces the surfactant density at 
which the micelles form; the critical micellar concentration (cmc) is reduced. Whilst if many micelles adsorb onto 
a single polymer the repulsive micelle-micelle interactions cause the polymer chain to swell; the radius of gyration 
may increase by a factor of two. Below, we propose a simple minimal model for necklaces of micelles adsorbed on a 
polymer chain. We then go on to calculate the free energy of such a necklace and use this to determine when they 
form and to look at the competition between free micelles and micelles bound to polymer. 

A polymer chain is a linear object, its topological dimension is 1 Q. So, two micelles adsorbed onto a chain cannot 
pass each other without one of them detaching from the polymer chain. Therefore the necklace of micelles along a 
polymer chain behaves as a 1-dimensional system. In the limit that the polymer molecule is very long the necklace is 
a 1-d bulk fluid. We assume that the micelle-micelle interactions are repulsive and so no two micelles can approach 
each other closely. This restriction on the relative positions of the micelles of the necklace restricts the configurations 
the polymer can adopt and so reduces its entropy. At high density of micelles the micelles are closely spaced along 
the polymer, that is the pieces of polymer between the micelles are quite short, and so these short pieces of polymer 
have to stretch in order to bridge the gaps between the micelles. This loss of entropy eventually limits the number 
of micelles which adsorb on a polymer and when the polymer chains are saturated with micelles, free micelles are 
formed. 

In the next section we derive from first principles a theory of objects with excluded volume interactions, which are 
adsorbed onto a linear polymer. This is the our main result, our results are not qualitatively new. They are however, 
semiquantitative, unlike previous work which was more qualitative. |^ First, the polymer-mediated interaction between 
two micelles is derived, yielding a potential of mean force for the intermicellar interaction. Then this is combined with 
a theory of 1-dimensional fluids to provide the equation of state of particles adsorbed onto a polymer. This is then 
used in section III to calculate the equation of state of a necklace. Combining the equation of state with a standard 
theory of micelle formation j^J^ then allows us to calculate how many micelles will form, both absorbed on polymer 
and free, as a function of the surfactant density. Finally, we discuss possible extensions of the theory and model. 

II. THEORY FOR NECKLACES OF MICELLES 

The system is a mixture of surfactant and polymer in a common solvent; the mixture is dilute, i.e., the fraction 
of the volume taken up by the surfactant and polymer is small. The polymer is composed of coils of size L, where 
L is proportional to the number of monomers of a polymer molecule [^|j9j. The coils are ideal, Gaussian, coils and 
if unperturbed the mean square distance between the two ends of a coil equals L. Thus L has dimensions of length 
squared, and L^^^ is the single relevant length scale for the polymer [^,^. We are assuming that the polymer's Kuhn 
length is much smaller than any length we consider explicitly. Our theory is entirely mesoscopic: the micelle is treated 
simply as a sphere and the polymer-micelle interaction is not considered explicitly but is treated phenomenologically 
via an association constant. This constant Afi is the difference in excess chemical potential between a micelle bound to 
a piece of a polymer chain and a free micelle. Here, excess means the chemical potential minus its ideal gas part: In p, 
where p is the (1- or 3-d) density. If Ap is negative the micelles will tend to adsorb on the chain This adsorption 
is due to part of the polymer chain lying near the surface of the micelle and interacting via an attractive interaction 
with it However, for simplicity we treat the polymer-micelle interaction as that between the micelle and a point 
on the polymer; we ignore the amount of polymer taken up in adsorbing to the micelle. With this simplification a 
micelle adsorbed onto a polymer coil is free to move along the entire length of the polymer L. As the polymer is a 
1-dimensional (1-d) chain of monomers the micelle acts as a particle on a wire. Even though the polymer is tracing 
out a an extremely convoluted path through 3-d space the configurational space available to a micelle adsorbed onto 
the polymer is only 1-d. 

If there is more than one micelle adsorbed on the same polymer coil then the micelles will interact with each other. 
We model the micelle-micelle interaction via a simple hard-sphere interaction: no two micelles' centres may be within 
a distance D of each other. The size of the micelles, D, is much less than that of a coil, L^/^, allowing many micelles to 
adsorb onto one coil. Above, we assumed that the polymer-micelle adsorption was pointlike, i.e., a negligible amount 
of the polymer is adsorbs on a micelle and we have just assumed that two micelles interact with each other out to 
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a distance D. For these two assumptions to be consistent the physical size of the micelles must be much smaller 
than the range at which two micelles interact. Essentially, we must have a micelle of diameter <C D which is quite 
highly charged so that no other micelle can approach within D . As micelles are typically 3-4 nm in diameter, 
D > 4nm. 

The limit of few micelles adsorbed onto a polymer chain, when micelle-micelle interactions are unimportant, is 
trivial (at our phenomcnological level of description). It is just the free energy of a polymer coil plus the excess free 
energy of a free micelle plus A/i plus the log of the 1-d density of micelles on the chain. The first of these is a reference 
free energy, which we need not specify, and the free energy of a free micelle and A/z are parameters of our model. 
However, at higher densities micelles on the same coil interact. They do so when the micelles are close together, i.e., 
when the pieces of polymer which connect adjacent micelles are short. The partition function of a piece of polymer 
chain with a micelle at each end is different from the partition function of the same length of chain without the 
micelles. In evaluating the partition function when the micelles are present, we must exclude those configurations in 
which the two micelles' hard spheres overlap. This is done in the following section and in section IIA we use that 
result to calculate the thermodynamic functions for a necklace of micelles in the limit that the coil is infinitely long. 
In IIB this result is used to determine free energy of a necklace at any density. 



A. The interaction between two micelles adsorbed onto the same polymer coil 

We consider two micelles adsorbed onto a coil and separated by a length of polymer L12. The probability of this 
separation of the micelles is proportional to the partition function of the chain at this separation. The two micelles 
have split the coil into 3 pieces, the one between them of length L12 and the two end pieces. The three parts are 
independent as the coil is ideal. Only the partition function of the centre piece is affected by the micelles; the end 
pieces are free to adopt any configuration but the the ends of the centre piece are each attached to a micelle and so, 
because of the micelle-micelle interaction, they must be at least D apart. The partition function of a chain of length 
L12 whose ends are constrained to be at least D apart is easy to calculate. The partition function for a chain with 
ends separated by r divided by its partition function when free, is just the probability distribution function (p for the 
ends of an unconstrained chain, times dr. Now, is |^ 

*(^-" = (2;;llj)"'»-(2£)- (" 

Thus, the partition function with the ends at least D apart, divided by the partition function for an unconstrained 
chain is simply 

An J (l,{Li2,rydr ^1- J exp{-x^)dx = exp[--w{Li2/ D^)]. (2) 

The seond equality defines w the effective interaction potential between two neighbouring micelles. The probability 
of the micelles being L12 apart, P{Li2/ D^), is directly related to the partition function of the chain when the micelles 
are L12 apart, 

J exp[-w[Li2/ D^)\dLi2 

for an infinitely long or ring polymer. We see that exp[— w{L 12 / D'^)] is acting as a Boltzmann weight and so w is an 
effective interaction; it is plotted in Fig. (|^). The interaction w is an athermal but soft repulsion. When two micelles 
are bound to the polymer less then L12 ^ D apart along a polymer chain then the chain has to stretch to cross the 
distance between the micelles, which must be > D. Conventionally, the Boltzmann weight is the exponential of minus 
an energy of interaction but here it is the exponential of minus the free energy of the piece of polymer coil. However, 
despite its unusual origin w behaves just as an interaction energy. The pair of micelles behave as a pair of particles 
on a wire which interact via a potential w. In the following section we use w to obtain the thermodynamic functions 
of a fluid of many micelles on one polymer coil. 



B. A 1-dimensional fluid of micelles on a polymer coil 



When many micelles are adsorbed onto the same coil they behave as a 1-d fluid. Formally we will take the 
L/D^ 00 limit, producing an infinite fluid. We do not consider effects due to the finite size of the polymer coils. 
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Each micelle's movement is restricted by the micelle in front and the micelle behind. As the system is effectively 1-d 
the micelles form a fluid at all densities; they never solidify or undergo a liquid-vapour transition Due to the 

convoluted path taken by the chain non-adjacent micelles can interact, indeed this interaction will cause the coil to 
swell and change from being Gaussian to being swollen |^Jl^. However, these interactions are weak in comparison to 
the interactions between nearest neighbour micelles and so will be neglected for the purpose of determining the free 
energy of the necklace of micelles. Each micelle always interacts with the micelles in front and behind it but collisions 
between different parts of the chain are relatively infrequent jl^ ]. 

Therefore we have a simple 1-d fluid. As the interaction between micelles is mediated via the polymer chain the 
interaction is strictly nearest neighbour. The free energy of the coil is simply a sum of the free energies of the pieces 
between the micelles and so is just a function of all the Li2's between neighbouring micelles. The statistical mechanics 
of 1-d fluids with nearest neighbour interactions is both well-studied and simple ||l3|. We start from the canonical 
partition function Zn{L) of N micelles on a coil of length L. As usual for 1-d systems, molecules fixed at each end 
of the system define its limits. So, with a micelle fixed at coordinate xq and one at a; at, Zn is 
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where the factor of 1/7V! is absent because we have put the micelles into a sequence: the i'th micelle is always the 
left-hand neighbour of the (i -|- l)'th micelle |13|. The units are defined so that the thermal energy kT is unity. 
Progress can be made if we move to the isothcrmal-isobaric ensemble [|l^,|l5[ ; introducing the pressure p which is the 
conjugate variable to the length L. The isothermal-isobaric partition function Zn(p) is 



ZNip) = J ZN{L)exp{-pL)dL. 



(5) 



Following the method of Takahashi we note that as the {N + l)'th particle defines the limit of our system 
xn+1 = L. Then if we define the new coordinates Vi = Xi — x^-i. 
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the partition function has reduced to a simple 1-d integration. The chemical potential p,id is, substituting N for 1 
in Eq. (|), 



P-id 



In / dr exp[— ti;(r) — pr], 

JO 



(7) 



which gives p^id as a function of the pressure p. We would like it as a function of the reduced 1-d density pid — ND"^ /L. 
The density pid is obtained from 



Pid 



\dpD^ 



(8) 



Using Eq. (||) in Eq. and then Eq. (||) yields p,id as a function of density pid; it is plotted in Fig. (^. Notice 
that the chemical potential is almost a linear function of density at high density. For future use we define an excess 
chemical potential pex = p-id — In pid. We now possess the chemical potential of a necklace of micelles as a function 
of the density of micelles along the polymer chain. In order to determine the density of micelles along a polymer as 
a function of surfactant density we require a model for the formation of micelles. 



III. MICELLE FORMATION 

Our model for micelle formation is quite standard 1^,01. We start by assuming that the surfactant monomers, the 
free micelles, and the micelles bound to polymer coils form an ideal ternary mixture. The only interactions considered 
are those between micelles bound to the same polymer; interactions involving free micelles or monomers are neglected 
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as the surfactant + polymer solution is dilute. The ternary mixture's density and composition are specified by three 
densities: the density of surfactant monomers pq = NqD^ /V, of free micelles pfm = NfmD^ /V, and of bound micelles 
Pbrn = NbmD^/V. Nq, Nfm and Nbm are the numbers of surfactant monomers, free micelles and micelles bound to 
polymer, respectively, and V is the volume. The 1-d density of micelles along the polymer chains pid is related to the 
density of bound micelles in the mixture pbm by 

Pbm = PpPld, (9) 

where pp = NpLD/V is a 'segment' density of the polymer. It is the number density of segments of polymer of size 
D^, times D^. Np is the number of polymer coils. Below, we treat the monomers and free micelles, and the micelles 
bound to polymer as constituting two different bulk phases. The polymer density pp just defines the relative volumes 
of these two parts of our system. 

At equilibrium the chemical potential of a surfactant is the same if it is a monomer or is part of either type of 
micelle. Then the chemical potential of a micelle composed of m surfactant molecules is m times that of a monomer. 
For the chemical potential of a monomer to be equal to that in a free micelle 

m(lnpo + Mo) = In (10) 
and for it to equal that in a bound micelle 

TO(lnpo + Mo) = Pldipid) + Mmic + = p,ld{pbm/Pp) + Prnic + A^, (11) 

where we usedEq. (|) to go from the middle to the right hand side expression. We have assumed that the number of 
surfactant molecules in a micelle is m for both free and bound micelles. Micelles are of course not monodisperse and 
so this is an approximation, although the numbers of micelles much larger or smaller than the most probable size is 
small I^J^. The constants po and pmic are the contributions to the chemical potentials of the monomers and micelles, 
respectively, of the interaction with the solvent. They are the change in free energy in excess of the Inp term, of the 
system when a monomer or a micelle is inserted . We have neglected the contribution from the momentum degrees 
of freedom. Thus, Pmic — ''npo is the difference between the chemical potential (minus the In p terms) of m surfactant 
molecules in a micelle and m surfactant monomers. 



The first equation for equilibrium between monomers and micelles, Eq. (10), is easily rewritten in the familiar form 



Pfm = Po" exp(m/lo - Prmc)- (12) 

The most convenient thermodynamic variable to work in is po which corres pon ds to controlling the chemical potential 
of the surfactant: Inpo + Mo- The second equation for equilibrium, Eq. (O), cannot be rewritten in the form pi,m 
equals a function of po because the the right hand side includes pid which is not a simple logarithmic or linear function 
of Phm- If Pp, pmic and Ap are specified, then Eq. ( |Tl| ) can be solved as a non- linear equation for pbm, at any value 
of Po. 

The total number density of surfactant is proportional to p = Po/'m + Pfm + Pbm- As all three of the densities po, 
Pfm and Pbm) are made dimensionless by multiplying by D^, the volume fraction of the surfactant is approximately 
equal to p times the physical volume of a micelle divided by D^. The physical volume of a micelle is the volume 
actually occupied by the surfactant molecules. As the physical volume is much less than the surfactant's volume 
fraction is much less than p. In Figs. (^) and (||) we plot the densities of monomers and micelles as a function of 
this total surfactant density. In Fig. (^ we have no polymer and there are only free micelles. The critical micellar 
concentration (cmc) is clearly around 10~^, although it is not that well defined as our aggregation number m is 
not that large, m = 20. In Fig. (^), the value of Pmic is unchanged from Fig. (^) but there is polymer present. As 
Ap is negative the polymer stabilises the micelles. This has two effects: micelles form at lower densities, the polymer 
encourages micelle formation, and micelles prefer to be adsorbed onto a polymer than to be free. So, at low densities 
where there are few micelles they are nearly all adsorbed onto polymer chains; the density of free micelles is very low 
below p ~ 3 X 10^^. As the density of micelles increases the density of micelles on the polymer chains becomes quite 
high and the interaction free energy of these micelles, p^x, is large. The polymer chains are stretching to accommodate 
the micelles and this stretching is reducing the configurational entropy of the polymer. This reduces the free energy 
reduction on adsorption and so we see that free micelles start to form. From Eqs. (p^ ) and ([III), it is easy to see that 
when /Xex = In Pp — Ap the numbers of free and bound micelles are the same. 

Notice that the density of surfactant at which free micelles appear is much higher in Fig. (^) than in Fig (^) because 
the monomer density is depleted by the formation of bound micelles. Free micelles appear at a po determined solely 
by Pmic, see Eq. ([L2|), but this corresponds to a higher p if there is competition for monomers. The po curve of Fig. 
(p|) looks as if there are two cmc's; it has two steps while in Fig. (0) it has only one. The first is where bound micelles 
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form and the rate of increase of po decreases sharply. Then as the density of bound micelles increases the polymer 
chains become crowded, micelle formation becomes less favourable, and the rate of increase of po increases. At yet 
higher density, free micelles form and the rate of increase of po decreases again. Of course, if the polymer density pp 
increases then the polymer chains saturate at higher surfactant densities and the second 'cmc' is pushed to higher 
values of p. 

The formation of a necklace of micelles causes the ideal polymer coil to expand for two reasons: it stretches the 
pieces of the polymer chain between the micelles and the micelle-micelle interaction results in the chain being swollen, 
like a polymer coil in a good solvent Without micelles the radius of gyration of the polymer Rq — 

At low densities of micelles, pid ^ 1 the radius of gyration is difficult to estimate as the micelles cover only a fraction 
of the chain. Therefore some parts, with micelles, repel each other but other, bare, parts do not. However, for pi^ > 1 
almost all the chain is covered with micelles. In order to minimise the stretching of the polymer the spacing between 
the micelles will be quite regular and near D. So, we can consider the chain as being composed of pidL/D^ segments 
of length ~ D. Then as the chain is swollen Rq ^ {puL/ D'^Y^^D. As D"^ this is much larger than L^/^. 



IV. CONCLUSION 



Perhaps the simplest possible model of polymer coils with necklaces of micelles has been proposed and studied. The 
simplicity of the model has allowed a detailed examination of its behaviour with only very minor approximations. This 
differs from previous theoretical work |^ which involved a highly approximate treatment of a much more detailed, 
and so complex, model. We have shown that the 1-d topology IJ] of a polymer chain implies that the necklace behaves 
as a 1-d fluid. This allowed us to determine the behaviour of the necklace using the highly developed theory for 1-d 
systems p^ . In combination with a standard theory for micelle formation it was then straightforward to calculate 
the densities of micelles formed, see Fig. (||). 

The model was kept simple for the sake of simplicity and clarity. For the purposes of comparison with experiment 
a more complex and accurate model is required. The theory presented here can be extended to deal with the most 
unrealistic features of the model. For example, the theory neglects the finite size of the polymer coils. Micelles are 
quite large, 4nm across, and so if the polymer coil is small there may be room for only a few micelles to adsorb 
without stretching the polymer so much that it is not favourable to adsorb another micelle. Then the fluid of micelles 
is not a bulk fluid. The finite size of the coil corresponds to a 1-d fluid between two walls, L apart. Confined 1-d 

|. The coils studied in experiment ||^ are in a good solvent and 
so are swollen. Therefore it would be worthwhile to extend our theory for ideal coils to self-avoiding coils. Although 
self-avoiding coils are more complex and difficult to deal with than ideal coils enough is known about their 

free energy when they are stretched to perform at least a scaling theory of a polymer with a micellar necklace. 
Experimental results pl] | on a semidilute solution of polymer clearly show the competition between polymer entropy 
and the free energy of adsorption which we have observed here. We hope that the theory developed here can be 
applied even to this system. 

It is a pleasure to acknowledge a careful reading of the manuscript by J. Doye . I would like to thank The Royal 
Society for the award of a fellowship and the FOM institute AMOLF for its hospitality. The work of the FOM 
Institute is part of the research program of FOM and is made possible by financial support from the Netherlands 
Organisation for Scientific Research (NWO). 



fluids have been considered extensively, see Ref. |16 
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FIG. 1. A schematic picture of a polymer coil with a necklace of micelles. The polymer chain is the black curve and the 
micelles adsorbed onto the chain are represented by dotted circles of diameter D. 

FIG. 2. The polymer mediated, effective interaction w between two micelles adsorbed onto a polymer coil, as a function of 
their separation L12. 



FIG. 3. The chemical potential of a 1-d fluid of micelles on a polymer coil, interacting via the potential w. 

FIG. 4. The density of surfactant monomers po (solid curve), and of free micelles pfm (dashed curve) as a function of total 
surfactant density p. mpo — Mmic = 45 and m = 20. There is no polymer present. 

FIG. 5. The density of surfactant monomers po (solid curve), of free micelles pfm (dashed curve), and of bound micelles 
Pfm (dot-dashed curve) as a function of total surfactant density p. m//o — Umic = 45, m = 20, Ap, = —15 and pp = 0.01. 
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